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Motivated by recent experiments with Josephson-j unction circuits we reconsider de- 
coherence effects in quantum two- level systems (TLS). On one hand, the experiments 
demonstrate the importance of 1/f noise, on the other hand, by operating at symmetry 
points one can suppress noise effects in linear order. We, therefore, analyze noise sources 
with a variety of power spectra, with linear or quadratic coupling, which are longitudinal 
or transverse relative to the eigenbasis of the unperturbed Hamiltonian. To evaluate the 
dephasing time for transverse 1/f noise second-order contributions have to be taken into 
account. Manipulations of the quantum state of the TLS define characteristic time scales. 
We discuss the consequences for relaxation and dephasing processes. 

1. Introduction 

The dynamics of quantum two-level systems has always been at the focus of interest, 
but recently has attracted increased attention because of the ideas of quantum computing. 
Several systems have been suggested as physical realizations of quantum bits allowing for 
the needed controlled manipulations, and for some of them first elementary steps have 
been demonstrated in experiments (see reviews in Ref. [p]]). A crucial requirement is the 
preservation of phase coherence in the presence of a noisy environment. In solid-state 
realizations, incl. those based on superconducting circuits, major noise sources are the 
electronic control circuit as well as material-specific fluctuations, e.g., in the substrate. In 
many cases one lacks a detailed microscopic description of the noise source, but frequently 
it is sufficient to model the environment by a bath of harmonic oscillators with frequency 
spectrum adjusted to reproduce the observed power spectrum. The resulting "spin-boson 
models" have been studied in the literature (see the reviews fHH), in particular the 
one with linear coupling and "Ohmic" spectrum. On the other hand, several recent 
experiments appear to be described by spin-boson models with more general couplings 
and different power spectra. In this article we, therefore, describe these extensions of the 
spin-boson model and analyze how more general noise sources and coupling schemes affect 
relaxation and dephasing processes. 

In the following Section we summarize what is known about relaxation and dephas- 
ing processes for spin-boson models with linear coupling. We add several extensions, 
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including a discussion how these processes are induced by specific quantum manipula- 
tions. We also present some exact results, shedding light on the important question of 
low-temperature dephasing. This analysis, as well as the rest of this article are appli- 
cable to any quantum two-state system subject to noise with appropriate coupling and 
power spectrum. However, since the extensions have been motivated by recent experi- 
ments with Josephson-junction qubits, we shortly present in Section 3 the physics of these 
systems. We indicate how, by proper choice of the operating point, one can tune to zero 
the linear coupling to the environment. Section 4 covers the extensions of the spin-boson 
model, which were motivated by the recent experiments. This includes noise sources with 
a variety of power spectra, with linear or quadratic coupling, which are longitudinal or 
transverse relative to the eigenbasis of the unperturbed Hamiltonian. For transverse 1/f 
noise we evaluate also second-order contributions, which turn out to be comparable to the 
lowest order. We conclude with a summary. The technical parts of our work, which are 
based on an analysis of the time evolution of the reduced density matrix of the quantum 
two-state system will be described in a separate publication. 

2. Spin-boson model 

In this section we review the theory and properties of the spin-boson model, which 
have been studied extensively before (see the reviews @,[|). A quantum two- level system 
coupled to an environment is modeled by a spin degree of freedom in a magnetic field 
coupled linearly to an oscillator bath with Hamiltonian 



H = Wctri + a z ^2c j (a j + a)) +TCh ■ (1) 
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The controlled part is 

ftctri = ~\ B z °* ~ \ B * °* = ~\ AE ( cos0 ° z + sin61 Gx ) ' ( 2 ) 
while the oscillator bath is described by 

Hh = ^ a ] a j ■ (3) 

3 

In the second form of Eq. (|J) we introduced the mixing ang le 9 = t&n~ l {B x /B z ), which 

depends on the direction of the magnetic field in the x-2-plane, and the energy split- 



ting between the eigenstates, AE = y-B^ + B 2 Z . In the standard spin-boson model it is 

assumed that the bath "force" operator X = J2j c j( a j + a ]) couples linearly to a z . 

In thermal equilibrium the Fourier transform of the symmetrized correlation function 
of this operator is given by 

S x {u) = ([X(t),X(t f )] + ) w = 2hJ(u)coth^ . (4) 
Here the bath spectral density has been introduced, defined by 

JH = l^S( U - Uj ). (5) 

n J 
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At low frequencies it typically has a power-law behavior. Of particular interest is the 
"Ohmic dissipation", corresponding to a spectrum which is linear at low frequencies up 
to some high-frequency cutoff ui c , 

71 

J (to) = — a Tiu®{(jj c — uo) . (6) 

The dimensionless parameter a reflects the strength of dissipation. In a physical system 
it depends on the amplitude of the noise and the coupling strength. Here we concentrate 
on weak damping, a < 1, since this limit is relevant for quantum-state engineering. Still 
we distinguish two regimes: the Hamiltonian-dominated regime, which is realized when 
AE ^> ak-gT, and the noise-dominated regime, which is realized, e.g., at degeneracy 
points where AE — > 0. 

In the Hamiltonian-dominated regime, AE ^> a k-gT, it is natural to describe the 
evolution of the system in the eigenbasis, |0) and |1), of H c t r i- 

|0> = cos 9 - |T) + sin 9 -\[) 

|1> = - sin ~ |T) + cos 6 -\i) . (7) 
Denoting by r x and t z the Pauli matrices in the eigenbasis, we have 
H = -\AE t z + (sin 9 r x + cos 9 r z ) X + H h . (8) 



2.1. Relaxation and dephasing 

Two different time scales describe the evolution in the spin-boson model The 
first is the dephasing time scale r v . It characterizes the decay of the off-diagonal elements 
of the qubit's reduced density matrix pit) in the preferred eigenbasis (0), or, equivalently 
of the expectation values of the operators r± = (1/2) (r^ ±ir y ). Dephasing processes lead 
to the following long-time dependence: 

(r±(f)> = tr [r ± p{t)\ « (t±(0)> e^ AEt ' h e"^ . (9) 

(Other time dependences will be discussed below). The second, the relaxation time scale 
Treiax; characterizes how the diagonal entries tend to their thermal equilibrium values: 

(r z (t)) - t z (oo) oc e ~* /Trclax , (10) 

where t z (oo) = tanh(Ai? /2feeT) . 

In Refs. m the dephasing and relaxation times for the spin boson model with Ohmic 
spectrum were evaluated in a path-integral technique. In the regime a 1 it is easier to 
employ the perturbative (diagrammatic) technique developed in Ref. f| and the standard 
Bloch-Redfield approximation. The rates are [] 

AE AE 1 

r re iax = = sin 2 9 —— coth — — — = sin 2 9 S x {oj = AE/%) , (11) 

ri zkrI ft 

2 k 1 1 

1% = r^ 1 = \ r rdax + na cos 2 9 = - r ro i ax + cos 2 9 S x {u = 0) . (12) 



1 Note that in the literature usually the evolution of (a z (t)) has been studied. To establish the connection 
to the results ( |ll|Jl^ ) one has to substitute Eqs. (|9|Jl0|) into the identity a z = cos 9 t z + sin 9 t x . Here we 
neglect renormalization effects, since they are weak for a<l. 
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We observe that only the "transverse" r^-component of the fluctuating field, proportional 
to sin (9, induces transitions between the eigenstates (0) of the unperturbed system. Thus 
the relaxation rate (11) is proportional to sin 2 9. The "longitudinal" r 2 -component of 



the fluctuating field, proportional to cos 9, does not induce relaxation processes. It does 
contribute, however, to dephasing since it leads to fluctuations of the eigenenergies and, 
thus, to a random relative phase between the two eigenstates. This is the origin of 
the contribution to the "pure" dephasing rate T* (0), which is proportional to cos 2 #. 
Relaxation and pure dephasing contribute to the rate 1^ = |r re i ax + T* . 

The last forms of Eqs. (|ll|) and (|12|) express the two rates in terms of the noise power 
spectrum at the relevant frequencies. These are the level spacing of the two-state system 
and zero frequency, respectively. The expressions apply in the weak-coupling limit for 
spectra which are regular at these frequencies (see Section ^ for 1/f noise). In general, 
fluctuations with frequencies in an interval of width r re i ax around AE/h and of width T v 
around zero are involved. 

The equilibration is due to two processes, excitation |0) — > |1) and relaxation |1) — > |0), 
with rates oc (X(t)X(l/)) u; =±AE/h- Both rates are related by a detailed balance 

condition, and the equilibrium value 7^(00) depends on both. On the other hand, r re i ax is 
determined by the sum of the two rates, i.e., the symmetrized noise power spectrum Sx- 

In the environment-dominated regime, AE <C a k^T, the coupling to the bath is the 
dominant part of the total Hamiltonian. Therefore, it is more convenient to discuss 
the problem in the eigenbasis of the observable a z to which the bath is coupled. The 
spin can tunnel incoherently between the two eigenstates of o z . One can again employ 
the perturbative analysis |J but use directly the Markov instead of the Bloch-Redfield 
approximation. The resulting rates are given by 

r re iax « B 2 J(27rhak B T) , (13) 
« l-nak^I '/h . 



In this regime the dephasing is much faster than the relaxation. In fact, we observe that as 
a function of temperature the dephasing and relaxation rates evolve in opposite directions. 
The a-dependence of the relaxation rate is an indication of the Zeno (watchdog) effect [[]]: 
the environment frequently "observes" the state of the spin, thus preventing it from 
tunneling. 

2.2. Sensitivity to the initial state 

To get further insight into dephasing phenomena we analyze some experimental scenar- 
ios. In particular, we discuss their sensitivity to details of the preparation of the initial 
state (cf., e.g., Ref. ||). Such a preparation is a necessary ingredient in an experiment 
probing dephasing processes directly. We observe that depending on the time scales of the 
preparation part of the bath oscillators (the fast ones) may follow the two-level system 
adiabatically, merely renormalizing its parameters, while others (the slow ones) lead to 
dephasing. We illustrate these results by considering the exactly solvable limit 9 = 0. 

Dephasing processes are contained in the time evolution of the quantity (r + (t)) obtained 
after tracing out the bath. This quantity can be evaluated analytically for 9 = (in which 
case r + = er + ), for an initial state described by a factorized density matrix pit = 0) = 
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Pspin <8> pbath (i-e., the TLS and the bath are disentangled). We find 

(a + (t)) = V(t) e- lAEt (a + (0)) , where V(t) = Tr(e^ >/V^ e 1 *^ 2 p bath ) . (14) 

Here the bath operator $ is defined as 

3 3 

and its time evolution is determined by the bare bath Hamiltonian, $(t) = e lHht §e~' lHht . 
To derive Eq. ([14]) we made use of a unitary transformation by U = exp (— icr z $/2). In 
the new basis the bath and the TLS decouple, which allows for the exact solution. 

The expression (|14"|) applies for any state of the bath (as long as it is factorized from the 
spin). In particular, we can assume that the spin was initially (for t < 0) kept in the state 
||) and the bath has relaxed to the thermal equilibrium distribution for this spin value: 
Pbath = Pt = Z7 1 e~ /3Wt , where = Hh + J2j c j( a j + a \)- 111 this case we can rewrite the 
density matrix as pbath = e**/ 2 /^"**/ 2 , with the density matrix of the decoupled bath 
given by pb = Z^e - ^ 15 , and obtain 

V(t) = P(t) = Tr (e- i$w p b ) . (16) 

The latter expression (with Fourier transform P(E)) has been studied extensively in the 
literature flj9|-|T2|| . It can be expressed as P{t) = exp if (i), where 



A POO 

K(t) = — / du 

ixn Jo oo 
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coth f — — — ) (cos cut — 1) — i sin cut 
\2k B T 



(17) 



For an Ohmic bath (^) at nonzero temperature and t > h/2k^T, it reduces to 

„ \ Sxico = 0) 2k B T , . 

ReK (t )« -t = -TTa—?—t . 18 

h h 

Thus we reproduce Eq. (^) with the dephasing rate given by ( fT2| ) in the limit 8 = 0. 

This result is easy to understand at the Golden rule level. Pure dephasing processes are 
associated with transitions during which only the oscillators change their state, while the 
spin remains unchanged. Since no energy is exchanged only oscillators with frequencies 
near zero contribute. An analysis of the Golden rule shows that the range is the inverse 
of the relevant time scale, in the present case given by T^. Since for an Ohmic model the 
spectral density J(u) vanishes linearly as u — > 0, exponential dephasing is found only at 
nonzero temperature T > 0. The situation is different in the presence of sub-Ohmic noise 
where, even at T = 0, the low-frequency oscillators may cause exponential dephasing (see 



the discussion in Sec. PO) 



For vanishing bath temperature, T = 0, one still finds a decay of (cr+(£)), governed by 
Re K (t) ~ —2 a \n(u c t), which implies a power- law decay 

(a + (t)) = (u c tr 2a e-* AEt / n (a + (0)) for t > l/u c . (19) 



This result is beyond the Golden rule. In fact all oscillators up to the high-frequency cutoff 
oj c contribute. While the power-law decay is a weak effect, it highlights the influence of 
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the initial state preparation and coherent manipulations on the dephasing. This will be 
particularly important in the sub-Ohmic and the super-Ohmic regime, as we will describe 
later. 

So far we discussed various initial states without specifying how they are prepared. We 
now consider some possibilities. A (theoretical) one is to keep the spin and bath decoupled 
until t — 0. In equilibrium the initial state which then enters Eq. (O) is the product state 
Pbath = Pb oc e~^ Uh . A more realistic possibility is to keep the bath coupled to the spin, 
while forcing the latter, e.g., by a strong external field, to be in a fixed state, say |f). 
Then, at t = 0, a sudden pulse of the external field is applied to change the spin state to a 
superposition, e.g., ^(IT) + II))- If the bath has no time to respond the resulting state at 
T = is |i) = ^(|T) + II)) <%> |gt)> where |gf) is the ground state of Ti^. Both components 
of this initial wave function now evolve in time according to the Hamiltonian (p. The 
first, |t) ® |gf), which is an eigenstate of ([[]), acquires only a trivial phase factor. The 
time evolution of the second component is more involved. Up to a phase factor it is given 
by ||) <S> exp(— i7i[t/%) |gj) where 7i[ = Tib — J2j c j( a j + a })- As the state |g|) is not an 
eigenstate of 7ij, entanglement between the spin and the bath develops, and the coherence 
between the states of the spin is reduced by the factor | (g|| exp(— z7ijt//i) |gj) | < 1. 

In a real experiment of the type discussed the preparation pulse takes a finite time, 
during which the oscillators partially adjust to the changing spin state. For instance, 
the (7r/2) a ,-pulse which transforms the state ||) — > ^(IT) + II))? can be accomplished by 
putting B z = and B x = A for a time span nh/2A. In this case the oscillators with (high) 
frequencies, huj 3> A, follow the spin adiabatically. In contrast, the oscillators with low 
frequency, hujj -C A, do not change their state. Assuming that the oscillators can be 
split into these two groups, we see that just after the (7r/2) x -pulse the state of the system 
is 4g (J|) <S> g|^> + ||) <8> gj 1 ^) <S> g 1 ^ where the superscripts 'h' and T refer to the high- 
and low-frequency oscillators, respectively. Thus, we arrive at a factorized initial state of 
the type as discussed above. However, only the low-frequency oscillators are factorized 
from the spin and contribute to dephasing during the subsequent evolution, while the 
high-frequency modes give rise to renormalization effects. Thus we reproduce the decay 
of (a + (t)) described above (i.e., for an Ohmic bath a power-law decay), however, with a 
preparation-dependent cutoff frequency %uj c = A. We may add that similar considerations 
apply for the analysis of the measurement process. 

2.3. Response functions 

In the limit 9 = we can also calculate exactly the linear response of r x = a x to a weak 
magnetic field in the x-direction, B x {t): 



xUt) = T &(t)(r x (t)r x (0) - r x (0)r x (t)} 
n 



(20) 



Using the equilibrium density matrix 



.r<i 



1 



1 + e~^ E 



IT) (tl®P t + e^ A£ ||) (||® Pi 



(21) 
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where pi oc exp(— ^H\) is the bath density matrix adjusted to the spin state ||) and 
similar for p±, we obtain the susceptibility 



xUt) = I- [P(t)e~ iAEt 



+ e-P AE P(t)e iAEt 



(22) 



The imaginary part of its Fourier transform is 
1 



2(1 + e- 



P(huo - AE) + e~ pAE P{hu + AE) 



-u 



(23) 



At T = and positive values of u we use the expression for P(E) from Ref. |T2j to obtain 



X, 



-P(huo - AE) = Q(hu - AE)- 



-270 



-2a 



2r(2a) 



- (hu) — AE) 



2a- 1 



(24) 



We observe that the dissipative part x'Lx has a gap AE, which corresponds to the minimal 
energy needed to flip the spin, and a power-law behavior as uj approaches the threshold. 
This behavior of x'x X (uj) parallels the orthogonality catastrophe scenario ||13|| . It implies 



that the ground state of the oscillator bath for different spin states, |gj) and |gj_), are 
macro scopically orthogonal. In particular, for an Ohmic bath (|6]) we recover the behavior 
typical for the problem of the X-ray absorption in metals [ |13| . 

The case of perpendicular noise, 9 = tt/2, has been treated in the literature (see e.g. 
Refs. p,14}|). Then, in the coherent (Hamiltonian-dominated) regime, Xx X ( u ) shows a 
Lorentzian peak around the frequency corresponding to the (renormalized) level splitting 
AE. The width of the peak is determined by the dephasing associated with the relaxation 
processes, i.e., by = r relax /2 as given in Eqs. (|T]]) and (|T2"D. 



3. Josephson junction qubits 

In this section we describe, by way of a specific example, some principles of quantum- 
state manipulations and the influence of various noise sources. For this purpose we con- 
sider a Josephson-j unction qubit based on a superconducting single-charge box displayed 
in Fig. la (for an extended discussion see Ref. |15|). It consists of a small superconduct- 
ing island coupled via Josephson junctions with effective coupling energy i^j($ x ) and the 
junction capacitance Cj) to a superconducting lead, and via a gate capacitor C g to a gate 
voltage source V g [TE\. The Hamiltonian of this system is 



n 



(Q-Qe 
2C 



£?j($ x ) cosy? with 



Q = - 



d 



i d(hcp/2e) 



(25) 



The excess Cooper-pair charge Q = 2ne on the island, which is an integer multiple of 
2e, and the phase difference ip across the Josephson junction are conjugate variables flT7 |, 
as indicated by the second part of Eq. (|25|) . The charging energy, with a characteristic 
scale Ec = e 2 /2C, depends on the total capacitance of the island, C = Cj + C g , and is 
controlled by the gate charge Q g = C g V g . We allow for a circuit with several junctions 
such that the Josephson coupling energy can be controlled by an applied flux <3> x . 




Figure 1. a) Josephson charge qubit; b) eigenenergies as functions of the gate voltage (for 
fixed <& x ) 



When the charging energy dominates, Eq Ej, it is convenient to rewrite Eq. (p5|) in 
the eigenbasis of the number n of Cooper-pair charges [IB 

W= £ (2ne - Q g ) 2 |w) (n| _ 1 ^ / |n) (n + ^ + |n + 1} (n| y (2g) 

In this limit it is obvious that eigenenergies form approximately parabolic bands as a 
function of the gate charge Q g , which can be labeled by n and are shifted relatively to 
each other by 2e (see Fig. lb). Near the degeneracy points, where two parabolas cross, 
the bands are split by an amount equal to Ej($ x ). (A similar picture had been proposed 
in a different context, that of a current-biased low-capacitance Josephson junction in 
Refs. ||H^P0|j .) If we concentrate on the vicinity of such a degeneracy point the system 
reduces to a quantum two-state system with Hamiltonian 

H = ~AE ch (V s ) a z - ^j($ x ) a x . (27) 

We note that by controlling the gate voltage we vary the difference in the charging energy 
of the two adjacent charge states, AE ch (V g ), and hence what corresponds to an effective 
magnetic field in the z-direction. Similarly the applied flux controls the Josephson cou- 
pling energy, corresponding to a control of the effective magnetic field in the x-direction. 
In several recent experiments |^I|-^3[1 quantum manipulations of Josephson charge qubits 
described by (pTD have been performed. 

Also for general relations between the charging energy and the Josephson coupling 
energy we have a well-defined quantum-mechanical problem. The system has a sequence 
of levels Ei(V g ,$x), which depend periodically on the gate voltage and the applied flux 
and can be controlled in this way. At low temperatures and for suitable applied ac-fields it 



may be sufficient to concentrate on the lowest two of these states [[24J] . The other extreme, 
Eq <C Ej, is realized in Josephson flux qubits, studied by other groups p5,|26]. Here the 
control has to be achieved solely by applied fluxes (possibly, different fluxes applied to 
different parts of the circuit). 
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The control of the qubit by applied voltages and fluxes is always accompanied by noise. 
As is clear from Eq. (p7|) a fluctuating control voltage 5V(t) couples to a z , while the flux 
noise £$(£) couples to o x . Since these noises are derived from linear circuits they are 
Gaussian and can be modeled by a bath of harmonic oscillators. They are completely 
characterized by their power spectra, which in turn depend on the impedance Z(u) and 
the temperature of the relevant circuits. 

For instance, the gate voltage fluctuations in the control circuit satisfy 

S v (u) = 2Re{Z t (uj)} hu coth [^^j ■ (28) 

They are governed by the total impedance of the circuit as seen from the qubit, Z^ x {uj) = 
iusCqb + Z~ x (uj\ which is determined by the external circuit Z(uS) and the capacitance 
C q b = (Cj + C g 1 )^ 1 of the qubit in the circuit. At the degeneracy point the voltage 
fluctuations couple linearly to the qubit, via Hi = — 2e(C g /C q b) 5V a z . In general, this 
noise is neither longitudinal nor transverse in the eigenbasis of 7i c tri, rather it is charac- 
terized by the angle 9y = tan~ 1 (£ , j/A£' c h). At low frequencies the circuit usually behaves 
as a resistor Z{iS) = Ry, and the results for Ohmic dissipation discussed in the previous 
section apply. They are characterized by the dimensionless parameter 

which is determined both by the strength of voltage fluctuations of the environment 
(oc Ry) and by the coupling of these fluctuations to the qubit (ex C g /C q b). We note 
that the circuit resistance is compared to the (typically much higher) quantum resistance 
Rk — h/e 2 ~ 25.8 kQ. At frequencies considered here the typical impedance of the 
voltage circuit is Ry ~ 50 Q. Furthermore, the effect of fluctuations is weakened if the 
gate capacitance is chosen small, C g <C Cj . Taking the ratio C g /Cj = 10~ 2 one reaches 
a dissipation as weak as ay ~ 10~ 6 , allowing, theoretically, for ~ 10 6 coherent single-bit 
manipulations. 

Similarly the fluctuations of the external flux can be expressed by the effective impedance 
of the current circuit which supplies the flux and the mutual inductance characterizing 
the coupling. Choosing the impedance purely real, Ri, we have 

R K (M dE 3 (Q x ) \* 

a * = TeT * a*. > ( 3 °) 



4R 7 \ $ o <9$ 



X 



and 6$ = ta,n~ 1 (AE c h/ Ej). For an estimate we take Rj ~ 100 Q of the order of the 
vacuum value. For M ps 0.01 — 0.1 nH and Ej/ks ~ 0.1 K we obtain a$ »s 10 -6 — 10~ 8 . 

In spite of these favorable estimates for the dephasing rates the actual experiments 
suffer from further, stronger noise sources. For instance, in charge qubits an important 



contribution comes from "background charge fluctuations". Their origin may be charge 
transfer processes between impurities in the substrate, which typically lead to a 1/ f power 
spectrum. This noise source can be modeled by further 2-state quantum systems p8| , |29| , 



however, when many of them couple weakly to the TLS they can be approximated again 
by an oscillator bath with appropriate spectrum. We will analyze this model in the next 
section. 
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There are ways to reduce the effect of the fluctuating environment |24], which can be 
demonstrated by considering the 2-state Hamiltonian (0): The idea is to tune the gate 
voltage to the degeneracy point V^o, where the difference in the charging energy vanishes, 
AE ch (V g0 ) = 0, and, furthermore, to tune the applied flux to a point where the effective 
Josephson coupling energy has an extremum, <9-Ej($ x )/<9$ x = 0. Hence the Hamiltonian 
is 

1 l dAE ch (V g ) l ggjCjj 2 f . 

H = --£j($ x0 ) a x - — 5V a z - - g ^ 2 5$ a x . (31) 

At this operating point the energy difference between the two relevant states depends 
only quadratically on the fluctuations in either the gate voltage SV(t) or the applied 
flux 8$>(t). On the other hand, one still can perform quantum-state manipulations by 
applying small-amplitude ac gate voltages, Kc(^), as is usual in NMR. (Only the noise 
at frequency AE/h couples through this channel; its effect is much weaker than that of 
the pure dephasing if 1/f noise is dominant.) By employing this idea the group in Saclay 
|24H has recently observed remarkably long dephasing times of the order of 1 //sec. This 
progress has been one of our motivations to study noise effects in more general situations. 

Another motivation to extend the analysis of the noise comes from the investigation 
of quantum measurement devices. Except when a measurement is being performed, the 
detector should be decoupled from the TLS as much as possible. Since it is hard to 
achieve complete decoupling, one usually turns to zero the linear coupling, but higher- 
order terms may still be present and affect the quantum dynamics. As an example we 
consider a single-electron transistor (SET), which can serve as a detector of the quantum 
state of a Josephson charge qubit [|T5| , |30| . The interaction between qubit and SET is given 
by 7imt = E- mt a z N, where N = J2i °i i°i ; 1S the number of electrons on the central island of 
the SET, and E- int denotes the capacitive coupling energy. The electron number N changes 
due to tunneling to the leads of the SET. Thus, in general the SET involves fermionic 
baths, but in some cases the problem can be mapped onto a spin-boson model: Consider 
the off-state of the SET with no transport voltage. In this case tunneling processes are 
suppressed at temperatures below the Coulomb gap of order Ec by the Coulomb blockade, 
and, classically, the charge of the central island is fixed at, say, N = 0. On the other 
hand, virtual higher-order processes (cotunneling) involve other charge states, e.g. N = 1, 
which makes N noisy. This situation can be mapped onto (the two lowest states of) a 
harmonic oscillator, with eigenfrequency equal to the Coulomb gap Ec, coupled to a 
dissipative bath. Since the interaction between the qubit and the oscillator depends on 
the occupation number operator of this oscillator, it corresponds to a nonlinear coupling 
in the formulation presented above. 

In this context we can also mention the recent proposal |H] of a SQUID-based device 



which allows switching off the linear coupling of a detector to a flux qubit, leaving also in 
this example only nonlinear coupling. 

These and further examples provide the motivation to study noise effects in more general 
situations than those covered by the usual spin-boson model. In the next Section we, 
therefore, investigate Ohmic circuit noise as well as 1/f noise, which couple in a linear 
or nonlinear way to the two-state quantum system. The noise can be longitudinal or 
transverse, which means in the eigenbasis of the TLS they couple to r z or r x , respectively. 



11 



4. Extensions of the spin-boson model 

The spin-boson model has been studied mostly for the specific case where a bath with 
Ohmic spectrum is coupled linearly to the spin degree of freedom. One reason is that 
linear damping, proportional to the velocity (Ohmic), is encountered frequently in real 
systems. Another is that suitable systems with Ohmic damping show a quantum phase 
transition at a critical strength of the dissipation, with o cr ~ 1. On the other hand, in 
the context of quantum-state engineering we should concentrate on systems with weak 
damping, but allow for general coupling and general spectra of the fluctuations. In these 
cases we find qualitatively novel results. For instance, a two-level system with sub-Ohmic 
damping (e.g., 1/f noise) still shows coherent oscillations, in spite of a general belief that 
it should localize in one of its states. Vise versa, in the super-Ohmic regime, where usually 
dephasing effects are believed negligible, we still find that manipulations with sharp pulses 
do influence dephasing. Finally we consider nonlinear coupling of the noise source to the 



qubit, which is important at symmetry points ||24j[ . Here, for definiteness we will consider 
only one source of fluctuations at a time. 

4.1. Sub-Ohmic and 1/f noise 

As mentioned above several experiments with Josephson circuits revealed in the low- 
frequency range the presence of 1/f noise. While the origin of this noise may be different 
in different circuits and requires further analysis, it appears that in several cases it derives 
from "background charge fluctuations" . In this case the noise of the gate charge may be 
quantified as Sq s (u) = ai/fe 2 /oj. Recent experiments |27] yield a>\/f ~ 1CT 7 -1CT 6 and 
indicate that the 1/f frequency dependence may extend up to high values, of the order 
of the level spacing of the TLS. 

One way to model the noise is to use an oscillator bath, where in the case of nonequi- 
librium 1/f noise the bath temperature should be treated as an adjustable parameter. 
For instance, a bath with J(u>) = (ir/2)ahuj s gives at uj <C k^T^/h the 1/f noise spectrum 

S x (u) = ^f. (32) 

with E\ij = 2nothuj s kBT\ ) . 

Such a model is a particular example of the so-called sub-Ohmic damping, defined by 
spectral densities 

J s ( u ) = (tt/2) haoj]~ s uj s with s < 1 . (33) 

Sub-Ohmic damping was considered earlier, for < s < 1, in the framework of the spin- 
boson model IHH , but did not attract much attention. It was argued that in the presence 
of sub-Ohmic dissipation coherence would be totally suppressed, transitions between the 
states of the two-level system would be incoherent and take place only at finite temper- 
atures. At zero temperature the system should be localized in one of the eigenstates of 
o z . The localization results from the fact that the bath renormalizes the off-diagonal part 
of the Hamiltonian B x to zero. However, this scenario, while correct for intermediate to 
strong damping, is not correct for weak damping. Indeed the NIBA approximation 0, 
which was designed to cover intermediate to strong damping, fails in the weak-coupling 
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limit for transverse noise. In contrast a more sophisticated renormalization procedure |32 
predicts coherent behavior. In the context of quantum-state engineering we are inter- 
ested in precisely this coherent sub-Ohmic regime. We will demonstrate that the simple 
criterion, which was used in Section |2| to distinguish between regimes of coherent and 
incoherent dynamics (i.e., the comparison of T* and AE), can be used also for sub-Ohmic 
environments. 

We consider the transient coherent dynamics of a spin coupled to a sub-Ohmic bath. If 
all the oscillators of the bath follow the spin adiabatically, the renormalized matrix element 
B x oc exp[— / du J '(lo) / 'u 2 ] is indeed suppressed, implying an incoherent dynamics. For a 
sub-Ohmic bath the integral diverges at low frequencies. However, a finite preparation 
time %/ A (cf. the discussion in Sec. 0) provides a low-frequency cutoff A for the oscillators 
that contribute to the renormalization, thus leading to a finite B x . The low-frequency 
oscillators, Tiujj -C A, contribute to the dephasing. Their effect is only weakly sensitive to 
the cutoff A, since the relevant integrals are dominated by very low frequencies. 

For longitudinal noise (9 = 0) with a 1/f spectrum one obtains for the function V(t) 
[Eq. (fUDI, with logarithmic accuracy p3p, RelnP(t) » -(E 1/f t) 2 \ ln(cu iT t)\/7fh 2 , where 
co>i r is the (intrinsic) infrared cutoff frequency for the 1/f noise.0 From this decay law one 
can deduce a dephasing rate, 

(34) 

Next we consider transverse noise, 9 = tt/2. As in the Ohmic regime we compare 
the energy splitting AE = B x to the pure dephasing rate T* (|34]) in order to determine 
whether the Hamiltonian- or noise-dominated regime is realized. In the former case, 
AE ^> for*, coherent oscillations are expected. This is confirmed by a perturbative 
analysis, similar to the one used to derive Eqs. (Pl , |i~2"l) , which gives in first order in a 








relax 


h AE ' 


r„ = 


2 1 relax 



(35) 

A more detailed analysis reveals that in the present problem higher-order contributions 
to the dissipative rates may be important as well. The reason is that, for 9 = 7r/2, the 
first-order relaxation rate is sensitive to the noise Sx(u) only at a (high) frequency AE/h 
where the 1/f noise is weak, whereas in second order the (diverging) low- frequency part 
of the spectrum becomes important. Employing the diagrammatic technique of Ref. M, 
we link the time evolution of the density matrix to a self-energy, E. We find that the 
Laplace transform of the second-order contribution E^(s) diverges as s — > iAE/h. Using 
first- and second-order terms in E we solve for the poles of t+(s) oc [s — iAE/h — E(s)] -1 
to describe the short-time behavior and find, with logarithmic accuracy 



E 2 1 E 2 

T v = a with a w — In - — ^— . (36) 

h AE 2vr huwAE 



2 In an experiment with averaging over repetitive measurements it may be determined by the time interval 
t av over which the averaging is performed In a spin-echo experiment the echo frequency may 

determine this frequency |27j]. 
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Thus, we see that the second-order correction to the self-energy changes the result for the 
dephasing rate considerably. As a consequence, for transverse noise the ratio T^/T^i^, 
which in first order takes the value 1/2, is shifted substantially towards higher values. It is 



interesting that in the recent experiments of Ref. [23] a ratio T^/T^i^ m 3 was reported. 
A more detailed description of the second order calculation will be presented elsewhere. 
Finally we note that Eq. ( Pop indeed confirms the assumption of underdamped coherent 
oscillations, as for re i ax , TiY^ <C AE in the Hamiltonian-dominated regime. 

Notice that for a sub-Ohmic bath with — 1 < s < 1 due to a high density of low- 
frequency oscillators the dephasing persists at Tb = 0: \(r + (t))\ oc exp[—a(u> s t) 1 ~ s /(s + l)] 
from which we read off the dephasing rate T* oc [a/(s + l)] 1 '^ 1_s - ) a; s (cf. Ref. p4f ). 

In the noise-dominated regime, AE <C for* , the dynamics is incoherent as we know 
from earlier work @,[|. Thus, our criterion for the coherent regime, AE ^> for*, is valid 
also for 1/f and sub-Ohmic environments. For < s < 1, the critical damping strength 



at which AE ~ for* coincides with that obtained in Ref. [32] and with the boundary of 
the applicability range of NIBA . 

4.2. Super-Ohmic case 

Next we discuss a super-Ohmic bath with, e.g., s = 2 and J s (u) = (tt/2) Tiauj 2 ju s . At 
zero temperature this leads to 

U s \ t J 

Thus, within a short time of order u~ l the exponent Re K(t) saturates at a finite value 
Re K(oo) = —2oiuj c /uj s . While usually this reduction of the off-diagonal elements of the 
density matrix is not denoted as dephasing (since coherent oscillations persist), from the 
point of view of quantum-state engineering it is a relevant loss of phase coherence, which 
introduces errors. We also note that this reduction strongly depends on the frequency 
cutoff uj c . Recalling that this cutoff is determined by the quantum-state manipulation 
(e.g., the width of a 7r/2-pulse) we observe that in the super-Ohmic regime the (weak) 
dephasing is strongly sensitive to the details of the preparation procedure. 

4.3. Nonlinear coupling to the noise source 

As motivated in the preceding section we study next a generalization of the spin-boson 
model to include nonlinear coupling of the noise source to the variables of the TLS. The 
Hamiltonian in the eigenbasis of the unperturbed system (with Pauli matrices denoted 
by r) thus reads 

1 X 2 

U = — AEt z + (sin 9 t x + cos 9 r z ) — + H h . (38) 

2 E Q 

Since the bath "force" operator X = J2j c j( a j + a ]) has dimensions of energy, for nonlinear 
coupling a new energy scale, Eq, appears. In specific systems this energy scale is a 
characteristic scale of the system. E.g., for a Josephson charge qubit coupled quadratically 
to the flux noise [Eq. ([H])] we have E = Ej. 

When considering a quadratic coupling we need to know the statistical properties of 
X 2 (t). While X is Gaussian-distributed, X 2 is not. However, for a first-order perturbative 
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analysis this fact is irrelevant and we can assume that X 2 is Gaussian as well, with a width 
that is fixed by the correlator S X 2 (u) = ({X 2 (t), X 2 (t')}) w . For baths with spectral 
densities which are regular at low frequencies it can be shown that this approximation is 
sufficient to describe the initial stage of the dephasing process for Y v t -C 1, and it provides 
a reliable estimate for the dephasing time. For singular spectra (e.g., for 1/ / noise) a more 
detailed analysis is needed, and non-Gaussian corrections may appear. Note that in order 
to evaluate the symmetrized correlator Sx 2 {u) in general one needs to know the full 
correlator (X(t)X(t')) and not only the symmetrized one Sx(w). For an Ohmic spectrum 
(H) the symmetrized and antisymmetrized spectra are related by detailed balance, and we 
find 

S X 2(uj) ~ a 2 hmax{(hu) 3 , (^bT) 3 } . (39) 



For 1/ f noise ( |32"D (assuming that the antisymmetric part diverges weaker than 1/u and 
can be neglected) we find 

S x2 (u) = -^ln-. (40) 

71 U U{ r 

Consider first the longitudinal coupling (8 = 0). For an Ohmic bath at low frequencies 
we have Sx 2 {u = 0) ~ a 2 H(k^T) 3 . Therefore, the pure dephasing rate is given by 

S X 2 (u = 0) a 2 {k B Tf 
* h 2 E 2 HE 2 ' 1 } 

For 1/ / noise, the correlator Sx* (u) also exhibits a 1/u divergence. Hence in analogy to 
the linear-coupling case the dephasing is governed by Re In V(t) = — [Ef,ft \n(u iv t) / 7ihE ] 2 , 
i.e., the characteristic dephasing rate is 

E 2 E 2 

r; = -^r\n^lL. (42) 
^ irri E§ u ir E 

In the transverse coupling case (6 = 7r/2) the relaxation and dephasing rates are again 
given by the noise power spectrum Sx 2 at frequency AE/h. In the lowest order we have 
Treiax = Sx 2 ( AE/K) / Eq for both Ohmic and 1/ / noise sources. In the case of 1/ f noise 
we should again take into account second-order corrections to the self-energy since they 
probe the divergent low-frequency part of the noise power spectrum. In this way we get 
T^ = b r re i ax with b ~ 1. 

4.4. Summary of results 

In the table we summarize our results for the relaxation and pure dephasing rates, r re i ax 
and r* of a two-state quantum system subject to the different sources of noise, with 
coupling which is longitudinal or transverse to the qubit (as compared to its eigenbasis), 
linear or quadratic, with an Ohmic (||), (|39l) or a 1/ f spectrum (|32|) , fl40"l). For the sake of 
brevity we put h = k-Q = 1 in the table and omit factors of order one in some expressions. 
Complete results can be found in the text. 
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linear co 


longitudinal ( ) 

pure dephasing) 
upling: 

H = ^AE Pz + X Pz 


transverse (±) 

(=>- relaxation + dephasing) 

W = ^AE Pz + X p x 


"Ohmic" 


r; = S x (uj = 0) ~ aT 


Trelax = Sx(AE) 
rrclax/2 


1/f 
quadratit 


(r + (t)> -e - ^//* 2 lnt 
r; = E 1/f \n l l 2 {E l/f /u iT ) 

2 coupling: 

H = \aE p z + -^X 2 Pz 


Trelax = Sx(AE) 

T ip = a r re i ax , a ~ 1 

w = \ae Pz + -^X 2 

Z tin 


"Ohmic" 


r; = ^ 2 (* = o)~a 2 ^ 


Trclax = T^^X^A-E) 

E 

Trelax/2 


1/f 




Trelax = S'x 2 ( A-E) 

T tp = b r rc i ax , 6 ~ 1 



5. Summary 

To summarize, motivated by recent experiments on Josephson-jucntion circuits, we 
have considered dephasing effects in a two-level system due to noise sources with various 
spectra (incl. Ohmic and 1/f), coupled in a linear or nonlinear, longitudinal or transverse 
way to the TLS. 

We have shown that the dephasing can be sensitive to the details of the initial-state 
preparation. For instance, a finite preparation time ~ K/A introduces an upper cutoff 
A on the frequency of environmental modes that contribute to dephasing, whereas the 
higher-frequency modes merely renormalize the parameters of the TLS. In particular, for 
an Ohmic environment we find a power-law dephasing at T = 0, sensitive to this cutoff 
frequency. 

We have also linked the mentioned renormalization effects to the behavior of a response 
funtion of the TLS, which exhibits features known for the orthogonality catastrophe, 
including a power-law divergence above a threshold. 

Noise with a 1/f spectrum we modeled by a sub-Ohmic bath and found a simple 
criterion for coherent behavior in sub-Ohmic environments. For transverse 1/f noise 
we found that, due to infrared divergences, one needs to take into account second-order 
contributions in a perturbative analysis. This can increase substantially the ratio between 
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the relaxation and dephasing times, which takes the value 1/2 in first-order calculations. 
These findings may contribute to the understanding of the experimental results of Ref. [pi], 
where a higher ratio (~ 3) was found. 

We have also analyzed the dynamics in the case of an environment with Ohmic or 1/ f 
noise power, which is coupled nonlinearly to the TLS. In the Ohmic case (e.g., for a single- 
electron transistor coupled to a Cooper-pair box) the dephasing rate scales as T 3 , and 
by biasing the SET at a special point it can be further suppressed to a T 5 -dependence. 
This should be contrasted, e.g., with a quantum point contact which, even in the off- 
state, couples linearly to the TLS and influences it strongly, Y v oc T. We reduced the 
description of a TLS coupled nonlinearly to a noise source with 1/ / spectrum to that of a 
TLS with linear coupling to an effective 1/f noise source. The results for the relaxation 



and dephasing times in this situation are relevant for recent experiments [24 
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